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. Abstract. The Sunyaev-Zel'dovich (SZ) effect of galaxy clusters is a tool to measure three quantities: Compton 

parameter, electron temperature, and cluster peculiar velocity. However, a major problem is non-removed contam- 
ination by astrophysical sources that emit in the SZ frequencies. This includes interstellar dust emission, infra-red 
(IR) galaxies, and radio sources in addition to primary Cosmic Microwave Background (CMB) anisotropies. The 
three former contaminations induce systematic shifts in the three SZ parameters. In this study, we carefully es- 
timated, both for a large beam experiment (namely Planck Surveyor) and a small beam experiment (ACT-like), 
■ the systematic errors that result if a fraction of the expected levels of emission from dust, IR galaxies, and radio 

sources remains non-removed. We found that the interstellar dust emission is not a major contaminant for the SZ 
, measurement. Unfortunately, the IR and radio source-induced systematic errors may be extremely large. In partic- 

^ ' ular the intra-cluster temperature and peculiar velocity will be determined inaccurately for Planck and ACT-like 

experiments, if only the frequency dependences are used for the cleaning. The Compton parameter is also affected 
. by the astrophysical contaminations. The systematic errors in this case were a factor of 2 to 5 times larger than 

' the expected statistical error-bar for Planck. For the ACT-like experiment, the statistical error-bars were larger 

than in the case of Planck by a factor of about 5, and therefore the systematic shifts remain within about 50% of 
the statistical errors. We have thus shown that the systematic errors due to contaminating astrophysical emissions 
can be significantly larger than the statistical errors, which implies that future SZ surveys aiming at measuring 
cluster temperatures and peculiar velocities will not be able to do so on their own without including additional 
• information like cluster shapes or follow-up observations. 
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' 1. Introduction these ex periments are interferoi netric arrays lik e 

AMIBA (iLo et al 200ll). AMI l)Kneissl et al. 200 1|) . 

The well-studied Sunyaev-Zerdovich (SZ) effect SZA ^ l|Mohr et al. 2002|) . Others are single dish 

rS ; llSunyaev fc Zeldovich 1971 ISunyaev fc Zel'dovich 1980) multi-frequency instrume nts like Plan ck surveyor^, 

^ . has been observed towards a few tens of known g alaxy SPT^, ACT^, OLIMPO (,Masi et al. 2003|) . SUZIE-II ^, 

■ ■ ■ ' clusters (see reviews by (jCarlstrom et al. 20021) and ACBAR^ MITO 

Birkinshaw 1999). It is a powerful tool for cosmological ^j^^ ^^^^^ ^1^^ b e used to characterise the 

and cluster studies, and when combmcd with other g^j^^^ ^j^^^g^^ themselves. llPointecouteau et al. 19981) 

observations (X-rays, optical, lensmg) the SZ effect proposed measuring the intra-cluster gas temperature 

allows us to measure cosmological parameters such as f^^^ relatiyistic corrections to the SZ. This is a partic- 

the Hubble constant and matter density in the uni- ^^^^^^ important issue in the context of future SZ surveys 

verse »^ (e.g. jMyers et al. 1997|), (|Gre^ et al. 200H), ^^-^^ X-ray counterparts will not be easily available 

^R^eese et al. 2002D , ( iBattisteUi et al. 2003D ). The SZ and since X-ray temperatures have not been measured for 

effect is independent of redshift and is therefore an clusters. This is because X-ray temperature determi- 

excellent tracer of large scale structure formation ^^^^^^ expensive, and one cannot expect that all the 
and evolution. This property has inspired groups 

to propose blind SZ surveys and investigate the i //astro.uchicago.edu/sza/ 

way they can probe the cosmological parameters 2 //www.rssd.esa.int/Planck/ 

e.g. (Bartlet t fc Silk 19941) , fBarbosa et al. 1996). 3 //astro.uchicago.edu/spt 
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Holder et al. 2000l). (Da Silva et al. 20001). 4 //www.hep.upenn.edu/~angelica/act/act.html 



^^^■^■^■^^^^^^^^^HaHK ' ; j j VV VV VV .lie: J^. UJ^ ^1111. ^UU^ dllg^lll^CLy' CLI^ U CLI^ lj . 11 bllll 

Kneissl et al. 200l|) . l)Xue fc Wu 200l|) '). As a conse- ^ //www.stanford.edu/~schurch/suzieJnstrument.html 



quence, several SZ experiments are either planned, ^ //cosmology.berkeley.edu/group/swlh/acbar/ 
under construction, or already observing. Some of ^ //oberon.romal. infn.it/mito/ 
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many clusters to be observed in future SZ surveys will 
have X-ray temperature determination. The SZ effect can 
also be used t o measure cluster radial pec uliar velocities as 
suggested by l|Sunvaev fc Zel'dovich 19800 . This provides 
a new distance free measurement of the cluster velocities, 
which is very interesting in view of future SZ surveys. 
However this needs measurement of the intra-cluster tem- 
perature. So far, it is only for a handful of known rich clus- 
ters that attempts have been made to extract the intra- 
cluster gas ten iperature i2) or to extract the peculiar ra- 
dial velocities llHolzapfel et al. 199'^ iLamarre et al. 1998t 
iBenson et al. 20031) . It is fair to sav that only upper limits 
have been obtained, and an actual temperature or peculiar 
velocity determination has still to be made. 

The SZ effect is a potentially powerful cosmological 
probe; however, it is not free of contamination. As a mat- 
ter of fact, SZ measurements are contaminated by other 
astrophysical sources that are mainly of two types: (i) due 
to our galaxy (free-free, synchrotron and dust emission 
from the milky way) or (ii) due to extra-galactic point 
sources (radio and infra-red (IR) galaxies), and even (iii) 
due to the Cosmic Microwave background (CMB) itself. 
The kinetic SZ and primary CMB fluctuations have the 
same frequency dependence, and therefore will a separa- 
tion of the CMB signal and the kinetic effect only be pos- 
sible through the relativistic corrections to the kinetic ef- 
fect. This would require both very high sensitivity (of the 
order O.l^K) and additional observing frequencies at the 
extrema of the relativistic corrections to the kinetic effect 
(near 200 and 500 GHz) and at its cross-over (near 300 
GHz) . The effect of the p rimary CMB anisotropies is well 
known (see for examp le l|Haehnelt &: Tegmark 1996|) and 
ijAghanim et al. 1997j) '). It induces a further uncertainty 
to the peculiar velocity and limits its accuracy. This effect 
intervenes mostly for large beam experiments like Planck, 
as small beam experiments are indeed less affected by the 
CMB which power is severely damped on a scale of about 
a few arcminutes. 

These contaminations may be monitored in the con- 
text of pointed SZ observations (i.e. towards known clus- 
ters) or follow-up observations of likely clusters. However, 
systematic follow-up observations of all the clusters in 
the SZ survey will be very time-consuming. Ignoring 
contaminations modifies the effectiveness of a survey by 
the loss of certain clusters and the appearance of some 
other "fake" artificial clusters. For future surveys accurate 
knowledge of the completeness and reliability is thus cru- 
cial for extracting cosmological information. The question 
of survey selection fu nction is being studied extensively 
for the Planck surv ey l|White 20031: iGeisbusch et al. 20041: 
ISchafer et al. 200^ : however, the systematic effect from 
the assum ed cluster structures still rema i ns to be studie d 
carefully ijBirkinshaw fc Lancaster 20041: iHansen 2004a|) . 
Contamination of SZ measurements by extra-galactic 
sources, especially radio sources, is n ot a new prob - 
lem (see for example reviews by l|Rephaeli 1995|) . 
Birkinshaw (1999)). It can be du e to radio em ission of 
the galaxies in the cluster itself l|Ledlow &: Owen 199& 



ICoorav et al. 1998t iLin et al. 2002|) or to foreground 
galaxies. More specificall v. the emissio n of radio sources 
can dilute the SZ signal. I Holder 2002() estimated the ex- 
pected dilution and its effect on the SZ power spectrum. 
IR dusty galaxies whose emission dominates at high fre- 
quencies may also contaminate SZ measurements. In par- 
ticular, gravitational lensing of dusty galaxies causes en- 
hancement of the confusion noise, which is lik ely to affect 
SZ observations of nearby clusters ijBlain 19981) . The effect 
of the point sources (both ra dio and dusty galaxies) hav e 
been recently revisited by ijWhite &: Maiumdar 200^ . 
who quantified it in terms of an equivalent noise and 
computed the as sociated power spectra. Also recently, 
l|Knox et al. 200^) investigated the effects of IR galaxy 
contamination on the statistical error-bars of the SZ pa- 
rameters. Contamination of SZ clusters by radio sources 
at low frequencies that can be monitored by the interfero- 
metric arrays remains potentially a very important source 
of errors for single dish experiments. On the other hand, 
multi-frequency observations with single dish experiments 
should help in solving the problem of contamination. 

In this study, we have focused on one category of SZ 
instruments planned for SZ surveying, namely the single- 
dish multi-frequency experiments. The SZ number counts 
that will be provided to us by these surveys will cer- 
tainly probe and constrain the cosmological parameters. 
However, such instruments are theoretically able not only 
to measure the SZ effect amplitude for each detected clus- 
ter (through the Compton parameter) but also its radial 
peculiar velocity and its gas temperature. We explored 
the capabilities of these future SZ surveys in terms of 
measuring the three cluster SZ parameters: Compton pa- 
rameter, intra-cluster gas temperature, and radial peculiar 
velocity in dependently of other observations. In a previ- 
ous work l|Aghanim et al. 20031) . we investigated the ef- 
fects of the cluster parameter degeneracies on the of the 
above-mentioned cluster SZ parameters in terms of error- 
bars. Here we were interested in the effects of the different 
sources of contamination on the measurement of these pa- 
rameters. The contamination did not affect the statistical 
error-bars strongly. 

The main question we are seeking to answer in this 
article is: How big are the systematic errors due to non- 
removable contamination? We show that the systematic 
errors are very important. They are by far the dominant 
source of trouble for the future single-dish SZ blind sur- 
veys when these surveys are used alone to estimate cluster 
parameters of unresolved clusters. We start by presenting 
our models for the major astrophysical contaminants in 
Sect. 2, and we then discuss the SZ parameter extraction 
technique and the sample of galaxy clusters used in our 
study in Sect. 3. In Sects. 4 and 5, we present the results 
for the three SZ experiments under consideration in our 
study, and discuss the results in Sect. 6. We finally offer 
our conclusions in Sect. 7. 
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2. Modeling the contaminants 

The contaminations strongly depend on the observing 
frequency. For galactic foregrounds the frequency depen- 
dance is given by the spectral energy distribution on the 
contaminants. For extra-galactic sources the frequency de- 
pendance results from the combination of the redshift evo- 
lution of the luminosity function with the spectral energy 
distribution of the sources. Moreover, the contamination 
levels, especially of the galactic foregrounds, also depend 
on the observed region of the sky and on the spatial dis- 
tribution of the signal (e.g. the power spectrum of dust 
emission) . Multi- frequency observations should help in re- 
ducing the contamination level through component sepa- 
ration, but we expect some contamination of the measure- 
ments to remain. This is rendered even more critical due 
to the rather large uncertainties in the emission models 
used to describe the contaminating sources. 

In order to address the question of statistical and sys- 
tematic errors on the cluster SZ parameter extraction due 
to non-removable contaminations, we computed the level 
of contamination at each source for a set of pre-defined fre- 
quencies. We based our study on an ensemble of planned 
SZ instruments which thus define our choice of frequen- 
cies. We distinguished between a large beam (~ 5 arcmin 
at best) experiment like Planck, and small beam 2 ar- 
cmin) experiments ACT-like. For the Planck instrument, 
we excluded from our analysis the 857 GHz channel which 
is totally dominated by IR emission from dust in the inter- 
stellar medium (ISM) and in high redshift galaxies. The 
ACT-like experiment had three observing frequencies 145, 
225, and 265 GHz. We assumed an optimistic sensitivity 
of 2 /iK for each channel and a beam of 1.7 arcmin at 145 
GHz, along with beam-size scaling as 1/frequency. 

2.1. Dust in the ISM 

Thermal emission from dust particles in the ISM is an im- 
portant source of radiation in the IR and submm domains 
{v > 100 GHz). It thus contaminates the measurement 
of the SZ effect at high frequencies. We therefore did not 
compute its contribution to the 30 and 44 GHz channels 
of Planck. Dust emission is usually modeled by a modified 
blackbody spectrum: 5^?""* = t x B^{Td) where S^(Td) is 
the blackbody spectrum, the dust temperature, and r 
is the emissivity. The emissivity is proportional where 
P is the spectral index. The values for (3 and Td are not 
unique but vary from dense clouds to diffuse gas. In the 
diffuse part of the sky (A^hi<5 10'^*' atoms/cm~^) the 
dust emission spectrum from the far-IR to the submil- 
limeter does not spa ti ally vary (Boulan ger et al. 199(i 
iFinkbeiner et al. 19991: lLagache~00 3') . Because of the 
difficulty of modeling dust emission properly, we instead 
com puted the observed dust emission spectrum as in 
l|Lagache 2003||) from a combination of COBE and WMAP 
measurements. We derived the mean spectrum of the dust 
emission correlated with HI in the diffuse part of the sky 
for |6| > 30°. This spectrum, which represents the far-IR 



to mm emission for iVHi=10^'' atoms/cm^, was used as a 
reference. The colors, i.e. the shape, of this spectruin are 
in very good agreement with l|Finkbeiner et al. 19991) . 

For the dust emission spatial distribution, our 
analysis relied on the power spe ctrum measured 
in HI jMiville-Deschenes et al. 20021) and in far-IR 
l|Gautier et al. 19921: iMiville-Deschenes et al. 2004j) . 
which is a steep power law power spectrum with no 
break. We thus used the same power spectrum oc k^^ 
that describes the dust distribution over a very large 
range of angular scales. 

To evaluate the fluctuation level of dust emission in the 
case of SZ experiments, we computed the dust emission 
power spectrum (in Jy^/sr) at each observing frequency 
using the following expression: 

P{k) = 4.8 10^ X {Nm/lO^y-^ X (fc/0.01)-3 X (Ba/Bioo)' 

where k is in arcmin^ ^, A^hi the HI column density in 
atoms/cm^ ^, and the -Ba/^ioo the color given at each 
observing wavelength A = hi>/c (in /im). In our calcu- 
lations, we used the column density associated with the 
cleanest 40% of the sky that corresponds to A'hi ^ 2.31 
10^" atoms/cm^. We finally computed the fluctuations of 
the dust emission in the beam of the considered experi- 
ments which are given by cr^ = J P{k)2-Kk dk Jy^/sr^. For 
a k~^ power spectrum, the previous integral was domi- 
nated by the lower limit of k, which is taken to be equal 
to 1/FWHM. 

The dust fluctuations cannot be removed and perfectly 
corrected for. In Tables ^ and |2] we show the level of con- 
tamination from dust fluctuations for Planck and for the 
ACT-like experiment for two cases: (i) where 30% of the 
fluctuations (i.e. 0.3(t) remain unremoved, and (ii) where 
100% remains. This 30% is the realistic level of contami- 
nation when independent tracers of dust emission (such as 
HI) are used to remove the dust contribution. The latter 
case (100%) is displayed to show how the maximum con- 
tribution affects the results. In both cases (30 and 100% 
fluctuations), the contamination is computed for the 40% 
cleanest region of the sky in terms of gas column density, 
as found outside the galactic plane and in bright molecular 
clouds. 

2.2. Infra-red or "dusty" galaxies 

Star-forming galaxies not only emit in the near-IR, 
optical, and UV domains, but also in the far-IR and 
submm. This emission is associated with the radiation 
that is absorbed and re-emitted by the dust in galaxies. 
For local galaxies, up to one third of the output emission 
lies in the far-IR domain. Star-forming galaxies can 
have much larger fractions of their emitted radiation 
in the far-IR (up to 90%). These galaxies are called 
"IR galaxies" in the following. The total emission from 
unresolved and/or faint IR galaxies is responsible for the 
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Frequency (GHz) 


30 


44 


70 


100 


143 


217 


353 


545 


30% Dust Fluctuations 
100% Dust lluctuations 




— 


3.96 
13.22 


13.09 
43.65 


33.76 
112.53 


144.69 
482.29 


704.14 
2347.15 


2218.54 
7395.14 


Poisson CIB fluctuations 
Poisson+Correlated CIB fluctuations 






311.27 
696.02 


583.41 
1304.54 


1445.88 
3233.09 


4426.90 
9898.85 


15279.8 
34166.7 


39082.8 
87391.8 


Radio source fluctuations 


1246.7 


1729.3 


2830.4 


3210.0 


3204.5 


3598.3 







Table 1. Level of contaminations (in Jy/sr) for Planck. The first two lines are for the dust emission fluctuations in the 
40% cleanest region of the sky, assuming a k^^ spectrum when 30% and 100% of the fluctuations remain unremoved. 
The second two lines are for the CIB fluctuations. We computed the case where the IR galaxies are Poisson distributed 
and the case where we also take the correlations between sources into account. In this case, a^^^ — v5>^^Poisson- 
The last line is for the fluctuations due to the unresolved radio, while resolved radio sources represent 12% to 15% of 
the background. 



Cosmic Infra-red Background (CIB) (Puget et al. 1996), 
which contains as much power as the optical and UV 
background s. The fluctuations of the CIB that are now 
detected llLaea che et al. 2000': 'Matsuhara e t al. 2000 l 
iMiville-Deschenes et al. 2002: Kiss et al. 20 01) are one of 
the major contaminants of the SZ measurements at high 
frequencies (> 100 GHz). As in the dust case, we will not 
compute the contribution at the 30 and 44 GHz channels 
of Planck. 

In order to evaluate the level of contamination by fluc- 
tuations of the CIB, i.e. unresolved dusty galaxies, we u sed 
the model of ijLagache et al. 20o3 iLagache et al. 20 04). 
which provided us with the number counts of IR galax- 
ies at each frequency of the SZ experiments (Planck and 
ACT-like). We flrst estimated the IR background emission 
from the dusty galaxies by integrating the number counts 
below a given detection limit Sj^. We then computed the 
rms fluctuations of the CIB due to unresolved sources in 
the beam of the SZ instruments. For this step, we con- 
sidered two cases, one where the sources follow a Poisson 
distribution and a second where t he correlation be tween 
sources is added to the shot noise. (ISong et al. 20031) have 
investigated the effect of correlations at high frequencies, 
and flnd it boosts the CIB shot noise fluctuations by a 
factor of 1.7 on average. Due to the lack of observational 
constraints or precise models for the correlations for all 
frequencies, we simply assumed that the correlated part 
of the CIB fluctuations is twice as large as the Poisson 
part for all frequencies and for all the experiments.^ It is 
worth noting that we did not take possible enhancement 
by lensing into account. Our estimate of the CIB contribu- 
tion is thus likely to be slightly underestimated especially 
for low redshift clusters. 

We estimated the contribution of the CIB fluctuations 
to the SZ measurements in two cases: one case where only 
1% of the dusty galaxies are resolved, which is t he case 
predicted for Planck surveyor l)Lagache et al. 2003|) . In the 
second case for the ACT-like experiment, we arbitrarily 
assumed 50% of the CIB is in resolved IR galaxies. The 



We have taken this factor 2 fr om the modeling of the clus- 
tering of the CIB anisotropies of iKnox et al. 20011). 



level of fluctuations are given for Planck and ACT-like in 
Tables Hand 121 respectively. 

2.3. Radio sources 

The radio emission of radio sources is dominated by 
synchrotron emission from relativistic electrons spiraling 
around magnetic fields in the sources. The synchrotron 
spectrum follows a power law v", with — 1 < a < 
—0.5. However, other types of sources exist such as those 
with "inverted spectra" (with a > 0). In general, the 
s pectra are not simple p ower laws with fixed indexes 
Herbig fc Readhead igg"^ . but have curved spectra in- 
stead; moreover some of the radio sources may be vari- 
able. It is therefore very hard to extrapolate their fluxes 
from one frequency to another. Additionally, radio sources 
are well ca talogued at low frequencies than ks to N VSS 
I Condon et al. 1998) and FIRST (White et al. 1997|) . but 
there are no all-sky surveys of extra-galactic sources at 
frequencies > 5 GHz apart from the WMAP observa- 
tions at 41 GHz jBcnne tt et al. 2003|) . As a result, prop- 
erties of the radio sources at high frequencies are poorly 
known. In the context of CMB observations, several stud- 
ies have aimed at predicting the radio source counts 
at frequencies higher than 5 GHz from k nown sources 
dToffolatti et al. 19981 ISokasian et al. 200 ih . These stud- 
ies allow us, in particular, to assess the level of contami- 
nation due to radio sources at the frequencies of interest 
for the CMB, where they are the dominant contamination 
between 30 and 217 GHz. However, one has to keep in 
mind that extrapolations of radio spectra from 5 GHz to 
high frequencies is quite a difficult task. 

In our study, we compute the level of contamination 
of unresolved radio sources per frequency by integrating 
the source counts below a flux limit, S1^^^°, determined 
for each frequency. Assuming Poisson distribution for the 
sources, we then evaluate the rms fluctuations due to un- 
resolved radio sources in the beam of each SZ experi- 
ment (Planck and ACT-like). We did not take the effect 
of the correlation of radio sources into account since it 
is rather weak. It is worth noting in the case of Planck, 
in particular, that the source fluctuations increase with 
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frequency although radio sources are getting fainter. This 
is due to the decrease of the beam-size with frequency 
from 33 arcmin at 30 GHz to 5.5 arcmin at 217 GHz. In 
practice, we use the radio-source counts that were kindly 
provided by L. Toffol atti and are based on the model of 
^Toffolatti et al. 199^ . This model was found to be in 
quite good agreement with recent observations of radio 
sources by WMAP at 41 GHz. It is thus likely to provide 
us with a good prediction of the radio counts at higher 
frequencies. As mentioned above, to estimate the level of 
contamination from unresolved sources we need the flux 
limit at each frequency for the de tection of radio sou rces. 
In the case of Planck frequencies, ijVielva et al. 200)^ pro- 
vided us with the flux limits 5'i'i^'° above which the radio 
sources are detected. They were computed using simu- 
lated observed sky taking all galactic and extra-galactic 
emissions into account, together with instrumental noise. 
These flux limits correspond to 5cr detections. As a conse- 
quence, the residual fluctuations computed in Tableware 
likely to be overestimated. 

Using the flux limits of ijVielva et al. 200)^ ranging 
from S-fi^f ° = 540 mJy (at 30 GHz) to S'i':^'° = 240 mJy 
(at 217 GHz), we found that about 12% to 15% of the ra- 
dio background was resolved into individual sources (de- 
tected at 5a) by Planck. The levels of contamination are 
summarised in Tabled For ACT-like, we investigated the 
effect of unresolved radio sources on the cluster param- 
eter measurement by assuming arbitrarily and optimisti- 
cally that 50% of the radio background is resolved into 
sources. The nur nbers are given in Table [3 Contrary to 
what is stated in I Knox et al. 200 J) , the radio source con- 
tamination always dominates over the fluctuations from 
IR galaxies up to 217 GHz when the latter are Poisson 
distributed, and up to 143 GHz when we take the correla- 
tions of IR galaxies into account. We did not compute the 
contributions at 353 and 545 GHz for Planck since they 
become quite subdominant. 



2.4. Synchrotron and Free-Free 

The synchrotron emission in our galaxy is due to the 
emission of free electrons spiraling around the lines of 
magnetic field. This emission is difficult to model be- 
cause of poor knowledge of the galactic magnetic field 
itself. Synchrotron emission i s expected to domin ate at 
low frequencies (< 80 GHz). ijBennett et al. 20031) argue 
that this emission in fact dominates over all other CMB 
foregrounds in the WMAP frequency bands. But at high 
galactic latitudes in the diff use sky (where the SZ obser- 
vations will be conducted), l|Lagache 2003|) show that the 
WMAP emission is dominated by the so-called "anoma- 
lous microwave emission" associated with the small inter- 
stellar transiently heated dust grains. This component was 
included in our measured dust spectrum. 

We did not include the free-free emission that is, 
in the diffuse sky at high latitudes, much smaller 
than the amplitude of the anomalous microwave compo- 



nent. We also ignore a p ossible contaminati on from ul- 
tra relativistic electrons ijEnsslin fc Kaiser^ OOO). which 
is at best marginallv detec table with future SZ sur- 
veys l Ensslin fc Hansen 20041) . 



3. SZ parameters 

The SZ effect is traditionally separated into two compo- 
nents according to the origin of the scattering of the same 
electrons 



AI{x) 



= A/thormal {x, y, Tg) -|- A/kinotic (x, T, Vp) (1) 

= y (gix) + St{x, Te)) - (3t h{x) , 



with X — hv/kTcrah and /q = 2(A:Tcinb)'^/('ic)^ where 
Tcmb = 2.725 K. The first term on the rhs of Eq. |(TJ is the 
thermal distortion ( Sunva ev & Zeldovich 1972*1 with the 
non-relativistic spectral shape 

gix)^ , (x'-^-a] , (2) 

and the magnitude is given by the Compton parameter 



y 



ctt 

me 



d/ Tip kTp 



(3) 



where rUe and rig are masses and number density of 
the electrons, ctt is the Thomson cross section, and the 
integral is along the line of sight. For non-relativistic 
electrons one has 5T{x^Tf.) = 0, but for hot clusters 
the relativistic electrons will slightly modify the ther- 
mal SZ effect llWright 1979|) . The spe ctral shape of 
5t(x, Tp) is easily calculated numericall v llRephaeli 199E ; 
DolgQv et a l. 2001 ; [S himon fc Rephaeli 2003 ; 

Sandoval- V illalbaz ofc Garcia-Colin 2003 
Itoh fc Nozawa 20041) . For the Planck frequencies 
thi s is easily impleme nted using the simple fits presented 
in l)Diego et al. 2003|) . 

The kinetic distortions ijSunvaev fc Zel'dovich 19*80(1 
in the last term of Eq. Q have the spectral shape 



h{x) = 



(e^ - ly 



(4) 



and the magnitude depends on = Up/c, the average line- 
of-sight streaming velocity of the thermal gas (positive if 
the gas is approaching the observer), and the Thomson 
optical depth 



T — I din 



(5) 



The minor thermal corrections to the kinetic effect are 
negligible with existing sensitivity. Thus, when the intra- 
cluster gas can be assumed isothermal one has y — 
TkTp/imeC^), while for non-isothermal clusters things are 
slightly more complicated (Hans en 2004) . 

Given the different spectral signatures of g{x),h{x), 
and 5T{x^Tf,)^ it is straightforward to separate the phys- 
ical variables and from sensitive multi-frequency 
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Frequency (GHz) 


145 


225 


265 


30% Dust Fluctuations 
lOOyo Dust fluctuations 


19.65 
89.17 


67.44 
224.79 


128.62 
428.74 


Poisson GIB Fluctuations 
Poisson+Gorrelated CIB Fluctuations 


1225.19 
2739.61 


5430.84 
12143.7 


10877.9 
24323.7 


Radio sources Fluctuations 


2265.42 


2873.19 


2821.53 



Table 2. Level of contaminations in Jy/sr for the ACT experiment. The first two lines are for the dust emission 
fluctuations in the 40% cleanest region of the sky when 30% and 100% of the fluctuations remain unremoved. The 
second two lines are for fluctuations of the 50% unresolved IR galaxies, when they are Poisson distributed and where 
we also take the correlations into account. In this case, (7^1^ = V^^^^'Poisson- ^^^^ represents the fluctuations 

due to unresolved radio sources when 50% of the background is resolved. 



observations. This separation is easily understood as fol- 
lows. To a good approximation the magnitude of the SZ 
effect gives the y parameter, which is measured at the 
maximum of g{x) (near 129 or 370 GHz). The velocity 
times optical depth, /3t, is measured at the zero-point 
of g{x) (near 217 GHz). Finally, electron temperature 
Te is measured using the zero point and maxima of 6t 
(near 190, 330, 475 GHz). Hence, under the assumption 
of isothermality, these observables can be combined to 
provide the physical parameters y,fp, and Te. It is thus 
clear that with at least 3 well-placed observing frequen- 
cies one can separate all 3 physical variables by using 
only the SZ observations. The first marginal determina- 
tion of all 3 SZ parameters was m ade for the h ot cluster 
A2163 (2) using data from SUZIE faolzapfel et al. 1997 ) 
and BIMA/O VRO llLaRoaue et al. 2002(l . However, as 
was shown in ("Aghanim et al. 2003 ?) additional observing 
frequencies are needed to break the intrinsic degeneracies 
between the SZ cluster parameters. 



3.1. The practical parameter extraction 

In order to find the SZ parameters for a given cluster 
we modified the publicly available SZ-parameter extrac- 
tion code SASZ. For a given "SZ observation" this code 
fast and accurately determines the central value for the 
3 SZ parai neters (y,Vr,, and Te), and corresponding la 



error-bars ( Hansen 20041) . These error-bars reflect solely 
the statistical error-bars due to the observational sensi- 
tivity (which for Planck or ACT-like is on the order of 
a few fiK for each channel). This parameter extraction 
code is based on the stochastic method of simulated an- 
nealing, and we chose the cooling scheme of SASZ in such 
a way that the central values of the 3 SZ parameters of all 
the clusters in our sample were determined to better than 
0.1% accuracy. 

3.2. The galaxy cluster sample 

To finally realistically answer the question of how big 
the systematic errors due to non-removable contamination 
are, one must first of all consider a large sample of clusters. 
We simulated clusters in bins of mass and redshift, using 



an extended Press-Schechter formalism that provides re- 
alistic temperatures, Compton p arameters, and peculi ar 
velocities (for details see Q) and ( Aghanim et al. 20031) ). 
In our study, we considered clusters with redshifts be- 
tween 0.2 and 3. Within this range of redshifts, not all 
clusters fill the beam of the instrument (i.e. 10i?c =beam- 
size, with Rc the cluster core radius), which sets a limit 
to the cluster mass range. We illustrate this point in the 
case of Planck, for which we set a conservative value of 
9 arcmin for the effective beam-size of the instrument for 
cluster observations. In this case, all clusters with total 
masses between 5. 10^^ and 10^^ solar masses and red- 
shifts z > 0.5 fill the beam. However, at lower redshifts 
z = 0.2(0.3 and 0.4) only clusters with masses above 
1.410^5(4. 10^^ and 8. lO^^) solar masses fulfill this con- 
dition. Specifically we consider 500 hot clusters with tem- 
peratures in the range 3.5keV < Tg < 7.5keV, central 
Compton parameter —4.7 < log{y) < —3.7, and peculiar 
velocity — 1600km/sec < Vp < 1600km/sec. For each clus- 
ter we thus have the true values, ytruc,7truc, and Utruc, 
from which one can calculate the true SZ signal, A/sz. 

Under the assumption that all contaminations have 
been removed completely, the total signal is equal to the 
SZ signal (ignoring the noise), A/totai = AJsz. Thereby, 
one can import the total signal to SASZ, which then finds 
the best fit parameters, j/dcrivcd: ^derived: and ^derived- 

For 

all the 500 clusters in our sample the central values of 
the parameters have been determined to better than 0.1% 
accuracy. 

Now, the interesting question is which central values 
are found for the SZ parameters when the contamina- 
tions have not been removed completely. For instance, 
if the dust has not been removed completely, then there 
will be an additional signal. A/dust- In that way the to- 
tal signal is A/totai = A/sz + A/dust- Believing this sig- 
nal is purely from SZ will induce a systematic error, 
I/derived — 2/true + ^ysystematic, but it wiU not significantly 
change the statistical error due to instrumental sensitiv- 
ity. We explicitly compared the statistical error-bars due 
to the instrumental sensitivity for slightly different cluster 
parameters, and found the effect on the statistical error- 
bars very moderate (a few percent). We will now address 
the question of how big these systematic errors are. We 
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show that these systematic errors for future planned SZ 
surveys can be significantly larger than the statistical er- 
rors due to instrumental sensitivity. 



4. Planck surveyor case 

As the first example we consider the all-sky Planck 
survey. The expected error-bars on the SZ parame- 
ters from Planck observations without taking the con- 
taminations into acc ount have already been studied 
l) Aghanim et al. 2003j) . Due to its wide frequency range, 
earlier predictions estimated that the Compton parame- 
ter should be easily measured and that the cluster pe- 
culiar velocity should be very marginally estimated with 
an error of about 500-700 km/s. In particular the impor- 
tance of non-Gaussian error-bars was discussed in detail in 
( Aghanim et al. 20031) . ^ n our study, we considered a sim- 
plification to the real error-bars. Instead of having differ- 
ent upper and lower error-bars, we simply took the largest 
of the two; that is, for all observables we consider TjJ^^"^ , 
with ST = (5Tup STiow = max((5rup, STiow)- 

In Fig. n] we present the expected Planck statisti- 
cal error-bars on y within this simplification. In the ide- 
alised case where no contamination is taken into account, 
Compton parameter y is found to an accuracy of 1% for 
the brightest clusters with y = lO"'^ '' and to 10% for the 
dimmer clusters in our sample with y =^ 10^^'^. The tem- 
perature, it is determined at best with a few keV, and 
more typically with 5 — 10 keV, statistical error-bars. The 
size of the statistical error-bars on Vp depends very much 
on the actual value of Vp, and goes from few hundred km/s, 
for the best case, to 2000-3000 km/s for the worst. It is 
worth noting that for our implementation of the extended 
Press-Schechter formalism (Sect. 3.2) we binned the data 
in cluster mass and redshift. This results in a quantiza- 
tion of the values (e.g. log(y)) and gives the impression of 
a binning in the figures. This has no impact on the results, 
since all we need is a sample of cluster parameters that 
roughly span the ranges expected in a standard cosmology. 

We next added the various non-removed contamina- 
tions to the signal to investigate their effects on the mea- 
sured parameters. More explicitly, with a pure SZ signal, 
A/inoasurcd = A/gj., One would derive the "true" SZ pa- 
rameters, whereas with a contaminated signal, A/obs — 
AJsz -|- A/contam One finds incorrect derived SZ parame- 
ters. Namely, the contamination causes systematic shifts. 
In order to quantify the importance of such systematic 
shifts, we compared it to the statistical error-bars, given 
above. We compared the value of e.g. Tdcrived — Ttmc to 
ST. Since the statistical error-bars from above are con- 
servative, then our findings represent lower limits to the 
magnitude of the systematic shift on the SZ parameters 
due to contamination. Our sample contains 500 clusters, 
and we considered (for Te, for example) the rms value. 
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over the sample, of the ratio Tdcrivcd — ?truc to ST, which 
is given by: 




We present the results in Figs. 11131 where we plot the 
central values of the derived parameters when adding the 
contaminations. In the case of interstellar dust on the 
40% cleanest part of the sky (filled circles), the systematic 
shifts are found to be quite small when 70% of the fluctu- 
ations are removed. The effect is small even when no dust 
fluctuations are removed. We found ey(dust) « ei,(dust) ~ 
0.15, and the temperature is only very marginally affected. 
We can therefore conclude that the contamination from 
the fluctuations of interstellar dust emission will not be a 
major problem for estimating the three SZ parameters {y, 
Te, and Vp) in the Planck case. 

The situation is completely different when we consid- 
ered the contamination due to the fluctuations of the IR 
background (CIB). In that case, we found ey(CIB) ft! 2.5, 
et,(CIB) 2, and found minor temperature shift. We also 
checked the effect of contamination by CIB Poisson fluc- 
tuations only. We found that the rms ratio between sys- 
tematic and statistical errors are about 0.5 for y and Vp 
and negligible for the temperature. 

The situation is even worse when we consider contam- 
ination from fluctuations due to the unresolved and unre- 
moved radio sources. In this case, we found a very large 
contamination of the Compton parameter ej, (radio) « 5, 
Et (radio) ^ (radio) ^ 1. As mentioned above, the level 
of radio fluctuations is likely to be overestimated since it 
was computed using flux limits for 5a detections. We re- 
did our analysis, assuming 20% of the radio sources to be 
resolved (less conservative than the 12-15% of Vielva et 
al. (2001)), in which case the Compton parameter are less 
affected by the radio contamination, Ej, (radio) » 3. When 
dust, CIB, and radio fluctuations are added there is a 
slight improvement of the determination of the Compton 
parameter. We flnd ey — 1.7; however, this does not im- 
prove determination of the temperature and velocities. We 
emphasize that this is an accidental partial cancellation, 
which cannot be counted on in general. 

The effects of the different contaminations on the SZ 
parameters are qualitatively easy to understand. They are 
illustrated with an ideal noise free experiment (0.1/iK) ob- 
serving at 100, 150, 217 and 270 GHz with about 1 arcmin 
resolution. First, a contamination from radio sources gives 
a positive contribution at low frequencies, where the domi- 
nating y parameter is determined, leading to a lower value 
for y (see triangles in Fig. ^ . Now, with this lower value 
for y, the high frequency signal seems too high, which can 
only be counterbalanced with a very low temperature (see 
triangles in Fig. [5)1. Finally, the velocity term reads as 
Vp/Tf,, so with a very low temperature, Vp is also forced 
to be very small. In reality, this leads approximately to 
Te ?» tip ft 0. 
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Fig. 1. The expected error-bars (red) on the Compton 
parameter for the Planck experiment, without contami- 
nation. The open triangles (blue) show the effect due to 
non-removed radio contamination. The open squares (pur- 
ple) are for correlated IR galaxies, and the filled circles 
(green) are when none of the dust emission fluctuations 
have been removed. 
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Fig. 2. Derived electron temperature in keV. Open tri- 
angles (blue) are for radio contamination, open squares 
(purple) are for Poisson-|-correlated IR galaxies, and filled 
circles (green) are for dust fluctuations. Statistical error- 
bars (of the order few to ten keV) are not shown in order 
to avoid confusion. 



Both IR galaxies and dust emission give a positive con- 
tribution at high frequencies. The latter can be compen- 
sated for in 3 ways, either 1) y parameter is larger (see 
squares in Fig. 0)), 2) temperature is lower (see squares in 
Fig- El J or 3) peculiar velocity Vp is slightly more negative. 
We see from Figs. I1I3I that the actual result is a rather 
complicated adjustment between these 3 possibilities. The 
lines of squares in Fig. 3 appear because of our quanti- 
zation in cluster masses and have no physical relevance. 
The systematic shift in peculiar velocity for a randomly 
observed cluster will be somewhere within the range of 
the symbols. The contaminating signal was added at each 
frequency for each cluster, rather than a population of 
contaminating signals; and therefore each cluster is repre- 
sented by only one symbol. We thus present an upper limit 



to the expected effect of contamination, and the effect for 
a given cluster can be smaller. 

In the Planck case, one also has to take the con- 
tamination from the CMB primary anisotropics them- 
selves into account. This contribution and its ef- 
fect were already studied and evaluated to an error 
of a few hundred km/sec ijHaehnelt fc Teemark 199^ 
^Agh anim et al. 1997t) . In particular, this was done in the 
context of a Wiener filtering analysis of the component 
separation taking a detailed simulation of the millimeter 
and submillimeter sky into account. The contamination by 
the primary CMB anisotropics slightly shifts the central 
values of the parameters without modifying the error-bars 
significantly. We therefore chose not to display it in the 
figures. 
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Fig. 3. Upper limit to the expected effect of contamina- 
tion in the derived peculiar velocity in km/s. Open trian- 
gles (blue) are for radio contamination, open squares (pur- 
ple) for Poisson-|-correlated IR galaxies, and filled circles 
(green) for dust fluctuation. Statistical error-bars (on the 
order of a few 100 to a few 1000 km/sec) are not shown, 
in order to avoid confusion. 



5. Future ACT-like SZ surveys 

A number of future experiments are planned to perform 
SZ surveys of a hundred to a few thousand square degrees. 
They expect to detect several thousand SZ clusters. We 
investigated the effects of the contaminations for those fu- 
ture ACT-like surveys that have a narrower frequency cov- 
erage but a smaller beam-size than the Planck surveyor. 

To be specific, we assumed that the future ACT-like 
experiment has three observing frequencies 145, 225, and 
265 GHz, with angular resolution of 1.7 arcmin at 145 
GHz (with beam-size scaling as 1 /frequency) , and sensi- 
tivity of 2 /iK for each channel. As a result of the smaller 
frequency coverage of such experiments, it will be more 
complicated for them to disentangle the SZ effect from 



Fig. 4. Compton parameter for an ideal experiment 
with 4 observing frequencies at 100, 150, 217, and 270 
GHz, and with 0.1/iK sensitivity. Open triangles (blue) 
for radio contamination; open squares (purple) are for 
Poisson-|-correlated IR galaxies. The very small statisti- 
cal error-bars without contamination are shown in red. 
Systematic shift from dust fluctuations is very small and 
is not shown as it would obscure the statistical error-bars. 

other sources which emit at the same frequencies. They 
will have a smaller beam-size of 1.5 arcmin as compared 
to the 5 arcmin Planck beam. This makes the contribution 
from interstellar dust smaller since it is dominated by large 
scale fluctuations. The CMB contamination will be quite 
negligible due to the sharp cut-off at small angular scales. 
On the other hand, those contributions due to fluctuations 
of CIB or unresolved radio sources will be larger than in 
the Planck case. However, more sources will be resolved; 
we thus set the fraction of resolved sources to an opti- 
mistic value of 50%. Moreover for ACT-like experiments, 
the number of clusters that fill the beam is much smaller 
than for Planck. Up to z = 0.4 all clusters with masses 
between 5. 10^^ and 10^^ solar masses will be resolved, i.e. 



4 



10 



Aghanim, Hansen & Lagache: Trouble for cluster parameter estimation from blind SZ surveys? 







>-2 



> 



^-4 



6 



-8 



"1 \ I r 



"1 r 



A 



A A 




A 

A A 



4 



5 

T 



true 



thousands km/sec. Therefore, even if a given contamina- 
tioii imphe ^ a complete failure in determining the central 
values, rtTSllrernains within the statistical error-bar. 

It is slightly different for the Compton parameter, for 
whichr expected er-for-barg ar 3 between 10 and 50%. Both 
duBt ^nd radio cOBtaminatioas will have a negligible im- 



pact (ej,(dust) ~ Ej, (radio) 



have a noticeable effect with systematic shifts in y of 

fand e„(CIB) = 0.4). It 
jal. 2004|) that including a 30 



up to 40% in the worst caies 
has been suggested ijKnox et i 



GHz channel (from complem sntary observations) will im- 
prove all these statistical error-bars (in particular for y), 
but as we saw above, thisj;hannel will be significantly 
:ontaminated by radio sources, which again makes y too 
inS'f'and pushes both and Vp towards zero). 
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Fig. 5. Derived electron temperature in keV for an ideal 
experiment. The very small statistical error-bars are 
shown in red. Open triangles (blue) are radio contami- 
nation, open squares (purple) are Poisson-|-correlated IR 
galaxies, and filled circles (green) are for dust fluctuations. 



exceed the beam size. The us e of adapted filterin g tech- 
niques like those proposed in l|Schafer et al. 200^ should 
help in reducing the contamination in such clusters. At 
higher redshifts up to z ^ 0.7, clusters with masses below 
2. 10^** fill the beam but they mostly represent the clus- 
ters close to the detection limit. It is only at high redshifts 
{z > 1) that more massive clusters (1. 10^^ solar masses) 
fulfill our condition on the beam. 

The few observing channels and moderate sensitivity 
of an ACT-like experiment makes determination of tem- 
perature and peculiar velocity very difficult. Statistical 
error-bars are typically much larger than the central val- 
ues themselves. For temperatures the smallest statistical 
error-bars are about 5 keV, and more typical error-bars are 
10-20 keV. For peculiar velocity the smallest error-bars are 
several hundred km/sec, and more typically it is several 
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Fig. 6. Compton parameter for ACT-like survey. The sta- 
tistical error-bars without contamination are shown in 
red. Open triangles (blue) are radio contamination, open 
squares (purple) Poisson-|-correlated IR galaxies, and filled 
circles (green) dust fiuctuations. 
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Fig. 7. Derived electron temperature in keV. Open trian- 
gles (blue) are radio contamination, open squares (pur- 
ple) are Poisson+corrclatcd IR galaxies, and filled circles 
(green) are dust fluctuations. Statistical error-bars (the 
smallest about 5 keV, but more typical error-bars are 
around 10-20 keV) are not shown in order to avoid confu- 
sion. 



Fig. 8. Upper limit to the expected effect of contamina- 
tion in the derived peculiar velocity in km/sec. Open trian- 
gles (blue) are radio contamination, open squares (purple) 
Poisson-|-correlated IR galaxies, and filled circles (green) 
dust fluctuations. Statistical error-bars (the smallest are 
about 500 km/sec, but more typical error-bars are 2000- 
6000 km/sec) are not shown, in order to avoid confusion. 



6. Requirements for SZ observations 

Our results for Planck show that the dust contamination 
on the 40% cleanest part of the sky is not a major source 
of error for the SZ parameter estimation even in the most 
pessimistic case when all the dust fluctuations remain un- 
removed. On the other hand, the fluctuating background 
from the IR or radio unresolved sources is an important 
issue. For the GIB fluctuations, we show that the Poisson 
fluctuations do not introduce major systematic shifts in 
contrast to the case where we take the correlations into 
account. The fluctuations from unresolved radio sources 
are the main sources of errors in the SZ parameters. This 
contamination dominates over the GIB fluctuations up to 



143 or 217 GHz when only Poisson or Poisson-t-correlated 
GIB fluctuations are present, respectively. 

In order to reduce the systematic effects on the SZ pa- 
rameters, we need to reduce the contaminations. A rather 
reasonable demand must be that for each contaminant the 
induced systematic error must be smaller than the statis- 
tical error, and hopefully much smaller. If we require that 
the ratio of systematic to statistical errors be smaller than 
1/2 for each contaminant, then we find that for Planck 
the radio contamination must be reduced by a factor of 
13, and the contamination from GIB must be reduced by 
a factor of about 5. For the radio sources, this is attain- 
able when ^ 60% of the radio background is resolved in 
individual sources, as compared to the 12-15% given by 
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the present detection limits for Planck. This corresponds 
to a flux limit of about 2 mJy at each frequency contam- 
inated by radio emission (from 3 mJy at 44 GHz to 1.5 
mJy at 217 GHz). For the GIB, reduction of the contami- 
nation by a factor 5 corresponds to resolving about 60% of 
the background with flux limits ranging between 0.67 and 
0.008 mJy from higher to lower frequencies. It is worth 
noting that due to confusion limit, at present only 50% of 
the GIB is resolve d in individual source s in blank fleld sur- 
veys at 350 GHz ( llfLa.sa.che et al. 20031) V This is achieved 
by the SGUBA instrument in small flelds (less than a few 
square degrees) and with a 10 arcsec resolution. 

In all cases (for Planck but also for AGT-like), the re- 
quirement to perform SZ parameter estimation from SZ 
observations is to observe on regions of the sky where the 
GIB and radio contamination levels can be pinned down. 
This is only possible by resorting to dedicated follow-up 
surveys, in radio and IR at higher resolution, to comple- 
ment the SZ observations. However, such follow-up ob- 
servations, in particular in radio, have to be done at fre- 
quencies close to the SZ measurements to avoid large er- 
rors from the extrapolation of spectral energy distribu- 
tions (Ricci et al. 2003). Both large flelds (a signiflcant 
fraction of the SZ surveys themselves) and angular resolu- 
tions (less than 3 arcsec typically) to beat confusion limits 
are mandatory. Future interferometers (e.g. ALMA) will 
have the required angular resolution but will not be able 
to cover large enough areas. More specifically, for ALMA 
to resolve 50% of the GIB in individual sources (number 
used for AGT-like experiments) at 230 GHz, the detec- 
tion limit is 0.1 mJy {5a). In such conditions, 138 days 
are needed to survey 1 deg^. It is even worse for resolv- 
ing 80% of the GIB; 96 days are needed for 10 arcmin^ 
with a detection limit of 0.02 mJy (5(t). For Planck, only 
the brightest sources need to be removed, which can be 
achieved with rather short observations. 

Gontamination by radio and IR point sources will be 
less of a problem for the interferometric SZ experiments, 
which naturally have better resolutions to monitor point 
sources. 

7. Discussion 

In our study, we have shown that the measurement of SZ 
parameters (y, Tg, and Vp), which is theoretically feasi- 
ble with multi- frequency observations of the clusters, is 
complicated in practice by the presence of unremoved or 
non-removable contaminations. In some cases, it is even 
made impossible. Even for very accurate observations of 
the SZ effect, i.e. with small error-bars, the contaminations 
induce systematic shifts in parameter values. These sys- 
tematic shifts can be extremely large, particularly in the 
case of two sources of contaminations: fiuctuations from 
unresolved radio sources and fluctuations of the cosmic in- 
frared back ground. Such cont aminations were recently ad- 
dressed by ijKnox et al. 2004|) . However, these authors es- 
timated the statistical errors on the SZ parameters rather 
than the induced shifts in the derived values. They were 



also studied by ijWhite fc Maiumdar 200^ . Rather than 
studying their effects on the SZ parameters, these authors 
estimated the level of noise introduced the IR and radio 
sources in SZ surveys, and compared it to the predicted 
SZ power spectrum. 

For the GIB and radio fluctuations, the systematic 
shift dominates the source of errors. Furthermore, the re- 
sults we obtained can become even worse, since we chose 
neither to take the amplification of th e GIB and radio 
emission by lensing effects into account ijBlain 1998fl . nor 
to include the effects of correlations for the radio sources 
and central radio source in clusters. Note that for AGT- 
like surveys we took an arbitrarily optimistic level of 50% 
for the unresolved IR and radio backgrounds, which can- 
not be achieved due to the resolution and the frequency 
of the experiment. 

The quantitative results we obtained for the system- 
atic shifts due to contamination by spurious sources of 
emission in the SZ frequencies are based on the hypothe- 
sis that the clusters under study fill the beam of the ex- 
periment. This will be the case for the major part of the 
Planck cluster catalogue; however, we do not expect this 
to be the case for AGT-like survey, for example, as the 
beam is much smaller 1.5 arcmin). Additional infor- 
mation provided by the spatial extension of the clusters 
helps in better constraining the optical depth r and thus 
the Gompton parameter itself. Such additional constraints 
will reduce the errors on the SZ parameter estimation for 
the next generation of experiments like AGT. They will 
not be as helpful in the case of Planck, where most of 
the clusters will fill the beam. Similarly, measurement of 
polarised emission from the clusters might help in better 
constraining r and thus y, provided the clusters do not fill 
the beam of the instrument. However, the polarised signal 
is weak and the measurement quite difficult. 

It seems from our results that the SZ observations 
alone will not provide us with the the three cluster pa- 
rameters. If we specifically want to estimate the peculiar 
velocities, we are forced to do it on known clusters. It 
is only on known objects that we might expect to prop- 
erly correct for the contaminations as was done e.g. in 
I Lamarre et al. on cluster A2163. However, such a 

procedure will consume time and effort. Another possi- 
bility might be to use extra information from scaling rela- 
tions and/or X-ray complementary observations. However, 
the use of X-ray observed scaling relations is unwise, 
partly because the weight from non-isothermal clusters is 
differ ent between X-r ay observations and SZ observations 
(e.g. l|Hansen 2004a(l ). and the use of numerical scaling 
relations is still far from being sufficiently reliable due to 
the complication in simulating gas dynamics. 

8. Conclusion 

The future blind SZ surveys will be contaminated by as- 
trophysical sources that can only be partly removed, such 
as interstellar dust emission, infra-red galaxies, and radio 
sources. Such non-removable contamination will induce a 



Aghanim, Hansen & Lagache: Trouble for cluster parameter estimation from blind SZ surveys? 



13 



systematic shift in the derived SZ parameters. We have 
shown that IR and radio source-induced systematic errors 
may be extremely large, virtually removing the possibility 
of measuring peculiar velocity and cluster temperature if 
using purely the SZ observations. Also the systematic shift 
in the Compton parameter will be significantly larger than 
the expected statistical error-bars. Therefore these con- 
taminations are potentially disastrous for future survey 
experiments, and must be considered very as very serious. 
On the other hand, the situation may be less desperate 
if the contamination levels are reduced by complemen- 
tary follow-up surveys. As mentioned in Sect.lHI the radio 
follow-ups should be planned at, or close to, the SZ fre- 
quencies and the radio source variability requires simul- 
taneous monitoring. GIB fluctuations need much better 
resolution than what is accessible now. One must there- 
fore be realistic about future SZ surveys, and consider 
them as excellent tools for identifying high redshift clus- 
ters and not as tools for measuring peculiar velocity or 
cluster temperature. 
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